
/* This is a bugy module */

#ifndef __SEPEWALD_H__
#define __SEPEWALD_H__

#include "complex.h"

#include "sepmisc.h"
#include "sepstrct.h"
#include "sepdef.h"

#define KAPPA .75

double sep_ewald_four(sepatom *atom, int kmax, double lbox, int npart);

double sep_ewald_short(sepatom *ptr, double cf2, double lbox, int npart, 
			  int ndim);
double sep_ewald_self(sepatom *atom, double lbox, int npart);

double sep_ewald(sepatom *atom, int kmax, double lbox, int npart);

void sep_3D_ewald_four(sepatom *atom, double kappa, int kmax, sep3D *sys);

void sep_3D_ewald_short(sepatom *ptr, double kappa, double cf, sep3D *sys);

void _sep_3D_ewald_short_bond(sepatom *ptr, double kappa, double cf, sep3D *sys);
 
void sep_3D_ewald_self(sepatom *atom, double kappa, sep3D *sys);

void sep_3D_ewald(sepatom *atom, double kappa, double cf_short, 
		  int kmax, sep3D *sys);


#endif
